Estimation of Coupling Constants of a Three-Spin Chain: 
Case Study of Hamiltonian Tomography with NMR 
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It has been shown that inter-spin interaction strengths in a spins-1/2 chain can be evaluated by 
accessing one of the edge spins only. We demonstrate this experimentally for the simplest case, a 
three-spin chain, with nuclear magnetic resonance (NMR) technique. The three spins in the chain 
interact through nearest-neighbor Ising interactions under site-dependent transverse fields. The 
employed molecule is an alanine containing three ^^C nuclei, each of which has spin-1/2. 



I. INTRODUCTION 

Fabricating a quantum system that would perfectly 
function as we desire is very challenging. Even with the 
most advanced nanotechnology it is still difficult to build 
a structure that would have the exact values of system 
parameters to make it work as we initially designed. This 
gives rise to the crucial necessity of system identification. 
In the context of quantum control, the system identifi- 
cation primarily refers to the identification of the system 
Hamiltonian, also known as Hamiltonian tomography. 

Yet, it is in general formidably hard to estimate the 
Hamiltonian: the number of necessary initial settings 
and measurements grows exponentially as the system size 
becomes larger. To make the problem of Hamiltonian 
estimation more feasible, various schemes to reduce the 
complexity and / or to minimize the effect of physical noise 
have recently been studied quite intensively. Examples 
are indirect, but efficient, schemes of Hamiltonian tomog- 
raphy of spin systems under limited access [H-O > and also 
an application of (classical) compressed sensing to the 
quantum setting 0] that greatly reduces the overall 
complexity. 

In the case of NMR, the interactions among spins in a 
molecule are usually determined by measuring all spins 
at once. But what if we are allowed to access only a 
single spin to reconstruct the whole Hamiltonian as in 
the above examples? In this paper, we report such an 
indirect Hamiltonian tomography of a three nuclei spin 
system with NMR as a minimum model. In the molecule, 
the spins effectively form a one-dimensional (ID) chain 
with nearest-neighbor interactions, which are of the Ising 
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type. We attempt to estimate the coupling constants by 
accessing solely the end spin pretending as if we had no 
knowledge about the interactions in advance. Then, we 
will make a comparison between the estimated coupling 
constants and the known values as estimated by standard 
methods. 

The analysis that we carry out is different from that 
presented in Q , where the coupling constants are es- 
timated from the energy eigenvalues obtained from the 
spectral peaks, which are obtained from long time evolu- 
tion of the spin at an end of the chain. In contrast, we 
obtain the coupling constants by fitting the time domain 
data without calculating the spectra in the present work. 
The materials presented in this paper therefore provide 
the first step toward the full verification of the scheme 
discussed in 0, @j. While data fitting is not computa- 
tionally efficient for large systems, we find that it is very 
suitable for the three-spin chain and more robust against 
relaxtion than [^, [^. Our finding thus paves the way 
to indirect estimation of a Hamiltonian of a small-scale 
noisy system where direct methods are not applicable. In 
such circumstances, it is impossible to obtain data over 
longer time periods, which is fundamental for obtaining 
sharp spectral peaks. 



II. THEORY 

In this section, we review how to estimate the spin- 
spin interaction strengths in a three-spin Ising chain with 
site-dependent transverse fields. The model in our mind 
is a three homonucleus molecule, such as alanine with 
three "'^^C nuclei. We use liquid-state NMR to control 
and measure the spins. 

The initial state is, thus, a thermal state 



Pth(r) = 



-Ho/ksT 



-Ho/kBT] 



(1) 



where 



Ho 
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with 1} = ^^I®!,!"^ =I®^®I,lf = I®I®^. 

Here T is the temperature, is the common Larmor 
frequency of the spins, and Ui is the zth component of 
the Pauh matrices. We drop the interaction terms among 
spins and chemical shifts of the spins temporarily since 
they are small enough compared with wq. We note that 
Eq. ([T]) is defined in the laboratory frame. 

A weakly coupled system develops according to the 
Hamiltonian 

n^ujiill+ UJ12II + ionll + J12II1I + J2zllll ■ (2) 

Here, uju and Jij characterize the transverse field of spin i 
and the coupling constant between spins i and j, respec- 
tively. We note that the Hamiltonian ^ is described in 
the rotating frames fixed to each spin. 
Now we evaluate the dynamics of spin 1 



Mk{t)^{il{t))^i:T[p{t)il], 

where k G {x, y, z\ and 

p{t) = (e-^«*)pth(r)(e-*«*)t 



(3) 



(4) 



is the density matrix of the system at time t. 

The dynamics of spin 1 without relaxations nor trans- 
verse field inhomogeneities is calculated and shown in 
Fig. [U when wii/(2T:) = 27 Hz for all i = 1,2,3. The 
coupling constants Ji2/(27r) = 53.8 Hz and J23/(27r) = 
34.8 Hz are taken from [9] as an example. 

It is clear that Fig. [T] is far from reality. We then 
include the effects of transverse relaxations via the op- 
erator sum representation e{p) = X)i=o ^jp^i where 
ELo ^1^^ = / [iS im. We take 

where 

Ao - -^p-*/T2(l) _u,-t/T2(2) _j_^-t/T2(3)^ 

A, 



Ao /, 

- A« {211) (* = 1,2,3), 



*/T2(l) _^g-t/T2(2) 

*/^^«) (z = 1,2,3), 



to represent the transverse relaxations. 

Finally, we take into account the inhomogeneity of the 
transverse fields 0Jn{x) as a function of position x. We 
assume that the inhomogeneity has a Gaussian distribu- 
tion 



P{uu{x)) 



1 



27rCT 



exp 



{uJu{x) - LOuf 



2a' 



(5) 



with variance a to be determined later. 

The coupling constants are estimated by comparing 
the spin dynamics obtained numerically with various 
pairs (Ji2,t^23) with the experimental data. Both the 
effects of the relaxations and the field inhomogeneity are 
taken into account in the numerical calculation. 
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FIG. 1: Ideal dynamics of the expectation values Mk{t) {k = 
x,y,z) of spin 1 when the initial state is a thermal state. 
Parameters uju/ {2n) = 27 Hz, Ji2/(27r) — 53.8 Hz, and 
J23/(27r) = 34.8 Hz are employed in Eq. 



III. EXPERIMENT 

We employ a linearly aligned three-spin molecule for 
demonstrating a Hamiltonian tomography through an 
edge spin. Our task is to determine the scalar coupling 
constants J12 between spins 1 and 2 and J23 between 
spins 2 and 3 by measuring only spin 1. 



A. Sample and Spectrometer 

We demonstrate a Hamiltonian tomography of a spin 
system with NMR. We employ a JEOL ECA-500 NMR 
spectrometer p^ . whose hydrogen Larmor frequency is 
approximately 500 MHz. We apply weak rf-fields to gen- 
erate transverse fields in the rotating frame of each spin. 

A 0.3 ml, 0.78 M sample of "C-labeled L-alanine (98% 
purity, Cambridge Isotope) solved in D2O, capsulated in 
a susceptibility matched NMR test tube [Ijl, is used. 
Three ^'^C atoms are linearly aligned in L-alanine. We 
label the carboxyl carbon spin 1, the a carbon spin 2, 
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and the methyl carbon spin 3. 

The scalar coupling constants are estimated from the 
spectrum obtained by Fourier transforming the free in- 
duction decay (FID) signal after a hard 7r/2-pulse is ap- 
plied for readout ^^j. Here, protons are decoupled using 
a standard heteronucleus decoupling technique (WALTZ- 
16) ^T^. The information extracted from the spectrum 
is summarized as follows. The Larmor frequency differ- 
ences are (a;o2— i^oi)/27r = 15.8 kHz and (woa— CLio2)/27r = 
4.4 kHz, where ujQi denotes the Larmor frequency of the 
spin i, for which the chemical shift is considered. Large 
differences in the Larmor frequencies compared with the 
scalar coupling constants justify the weak-coupling as- 
sumption made when the Hamiltonian ^ is introduced. 
The scalar coupling constant J13 between spins 1 and 3 
is on the order of 1 Hz [151] . which is much smaller com- 
pared to J12 and J23, and hence we can safely ignore it 
in our analysis. As a result, the Hamiltonian of alanine 
molecule is well approximated by Eq. 

Measured relaxation times are ri(l) — 15.5 s, Ti(2) = 
1.4 s, ri(3) = 0.9 s and T2(l) = 0.45 s, r2(2) = 0.23 s, 
72(3) = 0.63 s, where the argument labels the spin. The 
spin 2 has the shortest T2. In view of the fact that our 
data aquisition time to estimete the Hamiltonian is much 
shorter than any of Ti{i), we ignore the effect of Ti from 
now on. In contrast, we fully take the effect of T2{i) into 
account in our numerical calculations. 

The transverse fields are applied to spins by feeding 
oscillating currents with three different frequencies, cor- 
responding to the Larmor frequencies of the spins, to the 
coil. Their strengths are characterized by uju. 



B. Transverse Field Calibration 

We measure the dynamics of spin 1 in the presence of 
a;ii/(27r) — 27 Hz only, while other utu {i = 2,3) is set 
to 0, as shown in Fig. [21 The data was acquired in every 
0.004 s for < t < 0.6 s. The periodicity provides the 
information on the strength of the transverse fields, while 
the decay rate is determined by the relaxations and the 
field inhomogeneities. We find that the relaxations only 
are not enough to reproduce the decay, as demonstrated 
by the dashed line in Fig. [2j Both relaxations and field 
inhomogeneities must be considered to reproduce the de- 
cay rate. We obtain the variance a/utu = 0.05 in Eq. ([5]) 
by fitting the data. 



C. Results 

We measure the dynamics of spin 1 in two cases. 

In Case 1, the initial state is thermal and the trans- 
verse fields a;ii/(27r) = 27 Hz are applied to all the spins. 
The dynamics of the expectation values Mx,My, and Mz 
of spin 1 are shown in Fig. |3] as functions of time t. In 
Fig. [21 we see there are structures different from a simple 




FIG. 2: (Color online) Calibration of strength and inhomo- 
geneity of transverse field. Experimental results with a setting 
given in the text are shown. The dashed line shows the nu- 
merical result in which only the effect of relaxation is taken 
into account, while the solid line is the result in which both 
relaxations and inhomogeneities of transverse fields are con- 
sidered. 



sinusoidal oscillation which is expected without interac- 
tions. In other words, we obtain information concerning 
the interactions by measuring spin 1 only. It should be 
noted, however, that the spin dynamics is strongly af- 
fected by relaxations and transverse field inhomogeneities 
of all spins. 

In Case 2, we have chosen luh = 0, wi2/(27r) — 
a;i3/(27r) = 27 Hz. Since uiu = 0, we would not expect 
any dynamics in spin 1 if a thermal state were employed 
as an initial state. To avoid this problem, the initial state 
of spin 1 is prepared by applying a 7r/2-pulse along the 
y-axis to the thermal equilibrium state, while the initial 
states of spins 2 and 3 remain thermal. We again obtain 
the information on the interactions by measuring only 
spin 1 as shown in Fig. [31 



IV. ESTIMATION OF COUPLING CONSTANTS 

In this section, we pretend as if we do not know the 
coupling constants J12 and J23 and, instead, we estimate 
them by fitting numerically evaluated {Il{t)) with var- 
ious values of (Ji2,^23) to the experimental data. Al- 
though we defined the magnetization Mk{t) of the first 
spin as the expectation value (/^(f)), we temporarily as- 
sign Mk{t) to the experimental data while (/^(i)) to the 
corresponding numerical result to avoid confusion. Let 
us define the "distance" between the experimental data 
Mk{t) and the numerical result {ll{t)) by 



DkiJi2, J23) = J23)) - Mfc(t,)P- 

Here {tj} denotes the set of data acquisition points and 
{ll(tj, J12, ^23)) is the numerically evaluated expectation 
value of the fcth component of spin 1 at time tj with the 
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FIG. 3: (Color online) Dynamics of Mx, My, and Mz of spin 1 
are shown for Case 1, when all tJii/(27r) = 27 Hz and the ini- 
tial states of all spins are thermal. Experimental Results are 
shown by dots, and the solid lines are the numerical results, 
in which known values of coupling constants are employed. 



coupling constants (J12, t/23). In actual experiment, data 
was acquired in every 0.002 s for < i < 0.6 s. 

We do not make use of the spectra obtained by the 
Fourier transforms of time domain signals since the time 
window is not large enough to provide sharp peaks in 
the spectra. We can freely select the fitting window from 
t = to tw > ^0, where tg = 27r/Ji2 + 2tx j J2Z ^ 50 ms is 
the minimum time required for information to propagate 
from spin 1 to spin 3 through spin 2 and then propa- 
gate back to spin 1. Clearly there is an optimal value 
for t^, since too small provides too little data to be 
fitted, while too large makes relaxations and field in- 
homogeneities too significant. We estimate the coupling 
constants for three different values of and compare the 
results in the following. 



^ 0.8 



;d 0.6 
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FIG. 4: (Color online) Dynamics of Mx,My, and of 
spin 1 are shown for Case 2, when cjn = and tJi2/(27r) = 
cji3/(27r) = 27 Hz. The initial state of spin 1 is prepared by 
applying a 7r/2-pulse along the i/-axis to the thermal state, 
while those of spin 2 and 3 remain thermal. Experimental 
results are shown by dots, and the solid lines show the nu- 
merical results, in which known values of coupling constants 
are used. 



A. Case 1: cjn / 



In this case, the initial states of the three spins are pre- 
pared in thermal states and transverse fields u)u/{2'n) — 
27 Hz [i = 1, 2, 3) are applied to all three spins. 

3D and contour plots of the distances Dy{Ji2, J23) and 
Dz{Ji2, J23) with iw — 0.05 s are shown in Fig. O and 
Fig. [6l respectively. As expected, clear minima are found 
in these plots. 

We obtain the set (Ji2,^23) that mini- 
mizes -Dj,( J12, ^23) and Dz{Ji2,J2?,) for different 
tw = 0.05, 0.1, 0.2 s. They are summarized in TableHl 
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30 32 34 36 38 



J23/27r [Hz] 

FIG. 5: 3D and contour plot of Dy{Ji2, J23) for Case 1 when 
tvr = 0.05 s. A clear minimum can be seen. The distance 
between two neighboring contours in the contour plot is 0.01. 




28 30 3? 34 36 



./23/27r[Hz] 

FIG. 6: 3D and contour plot of -Dz(Ji2, J23) for Case 1 when 
tw ~ 0.05 s. A clear minimum can be seen. The distance 
between two neighboring contours in the contour plot is 0.01. 



TABLE I: Estimated coupling constants in Case 1 for various 
window size tw. 



Case 1 


tw [s] 


Ji2/27r [Hz] 


J23/27r [Hz] 


Known values 




53.8 


34.8 




0.05 


53 


35 


Dy 


0.1 


55.5 


36 




0.2 


55.5 


36.5 




0.05 


55 


33 




0.1 


54.5 


30.5 




0.2 


56.5 


25.5 



tance Dx{Ji2,J2j,) with different time windows tw = 
0.05,0.1,0.2 s. Table ITTl summarizes the results. 

TABLE II: Estimated coupling constants in Case 2 for differ- 
ent window 

Case 2 [s] Ji2/27r [Hz] J2s/2ti [Hz] 



Known values 




53.8 


34.8 




0.05 


56.5 


37.5 




0.1 


58 


38 




0.2 


59 


38.5 



B. Case 2: tJn = 

We take loh = and a;i2/(27r) — CL>i3/(27r) = 27 Hz in 
Case 2. To introduce nontrivial spin dynamics to spin 1, 
a 7r/2-pulse along the y-axis, Y = exp(— i7r/^/2), is ap- 
plied to spin \ &t t — Q after the thermal state has been 
prepared. 

Figure[7]shows the 3D and contour plots of the distance 
DxiJi2, J23) with = 0.05 s. As expected, a unique 
minimum can be found in the figure. 

We obtain the pair {Ji2,J23) that minimizes the dis- 



C. Estimation 

Regardless of the choice of or {D^, Dy, D^), the es- 
timated pair ( J12, J23) is consistent with each other both 
in Cases 1 and 2. It seems, however, that the smallest 
tw = 0.05 s yields the best results when we compare them 
with the known {Ji2,J23) obtained by different means. 
We have also confirmed that a smaller value, = 0.02 s, 
is not large enough to estimate the coupling constants 
reliably as shown in Fig. |51 where the profile has a sharp 
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.hill-K [Hz] 

FIG. 7: 3D and contour plots of Dx{Ji2, J23) for Case 2 with 
time window tw = 0.05 s. A clear minimum can be seen. The 
distance between two neighboring contours in the contour plot 
is 0.01. 



minium along the Ji2-direction but is almost flat along 
the J23-direction. This bahavior clearly shows the signif- 
icance of the time to ~ 0.05 s defined previously. The 
effect of J23 does not manifest itself yet in the behav- 
ior of spin 1 for a short time less than Iq. On the other 
hand, for ^ to, relaxation and field inhomogeneity are 
less serious yet, and the data produces an excellent re- 
sult, while the results provided by a larger suffer from 
these effects. 



V. SUMMARY 

We have successfully demonstrated for the first time 
that indirect Hamiltonian tomography is possible in 
NMR setup. As long as the system is small enough for 
efficient data fitting, the estimated values are surpris- 
ingly close to the real ones, given the substantial amount 
of noise and inhomogeneities in the system. This paves 
the way towards the identification of spins and couplings 
which are off-resonant or would usually the drowned 
by background noise. While the methods of |^ rely 
on Fourier analysis, which is only applicable in systems 




J23/27r[Hz] 



FIG. 8: 3D and contour plots of Dy{Ji2, J23) for Case 1 with 
time window t„ = 0.02 s. It shows that the time is not 
large enough to estimate the coupling constant J23. The dis- 
tance between two neighboring contours in the contour plot 
is 0.0028. 



clean enough for sufficiently long time data aquisition, 
our method can be applied in more noisy cases. 

We have shown that there is a competition in the ob- 
served evolution between amount of data acquired and 
the amount of noise coming in. It seems optimal to 
choose rather short data acquisition times in order to get 
a good agreement of the estimated couplings with their 
real values. 
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